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The advent of reasonably fast and inexpensive Very Large Scale Integration (VLSI) 
components offers an opportunity to reconsider the merits of composite-code-uplink 
ranging systems , abandoned years ago . It is shown in this article that a ranging receiver 
with a sufficient and reasonable number of correlators is competitive with the current 
sequential component ranging system and may outperform that system by some 7.5- 
2.5 dB. The optimum transmitter code, the optimum receiver, and a near-maximum- 
likelihood range-estimation algorithm are presented. 


I. Introduction and Background 

Prior to 1973, planetary and lunar spacecraft ranging sys- 
tems at JPL utilized a transmitted uplink code made by com- 
bining binary clock and pseudonoise sequences in a majority- 
vote logic [l]. 1 The “composite-code-uplink” ranging receiver 
consisted of one or two channels that correlated the trans- 
ponded signal with combinations of the clock and each sepa- 
rate component sequentially through each successive symbol- 
delay to determine the precise delay on each of the compo- 
nents. Since transmitted power was distributed among the 
various clock and pseudonoise code components but the 
receiver was sensitive to only one component at one phase at 
a time, acquisition time was longer by about a factor of 
16 than if the receiver could have processed all the received 
power during the acquisition time. 


1 R. W. Tappan and R. C. Tausworthe, “DSIF Technical Description, 
Planetary Ranging Equipment, Mariner-Mars 1971 Configuration,” 
JPL Document TD505943A (internal document), Jet Propulsion 
Laboratory, Pasadena, California, February 4, 1972. 


Consequently, once sufficient analysis and precautions had 
been taken to ensure that uplink ranging sidebands would not 
interfere with the spacecraft command system, a “sequential- 
component-uplink” ranging method [2] was devised that is 
still being used today. The term “sequential,” in this case, 
refers to the transmitted code, which is a time series of square 
waves of successively shorter wavelengths. The receiver, still a 
few -correlator device, is programmed to acquire the sequen- 
tially transmitted components one by one. The newer method 
had a 16:1 acquisition-time advantage over the older scheme 
because it could utilize all the transponded power for each 
component during the acquisition time. The necessity to pro- 
gram the uplink code and to program the receiver to switch 
components at the proper round-trip-light -time interval was 
a disadvantage compensated by the signal-to-noise advantage. 

During a transition period, both ranging techniques were 
used. The older, composite transmitted code method was 
referred to as the “r” system, and the newer, sequential com- 
ponent code method was called the “/z” system. These names 
were dubbed by Robertson Stevens, now the Chief Engineer 
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of the Deep Space Network, one day while strategizing at the 
blackboard in a design meeting, groping for a notation to dis- 
tinguish the two. They derive from the initials of the then- 
purveyors of the two systems: the author (guess which one) 
and Warren Martin. Usage of these designations has decayed 
over the years because the r system is no longer extant. The 
designations are reinstituted in this article for brevity in refer- 
encing the two schemes. 

The r planetary composite transmitter code was generated 
by combining a clock square wave with a majority-vote logic 
of 5 pseudonoise sequences in an exclusive -OR fashion. The 
components had symbol-periods of 2, 7, 11, 15, 19, and 23, 
for a total code period of N = 1 ,009,470. This, when clocked 
at a symbol-period t 0 of about 1 /isec, gave a repetition per- 
iod of about 1 sec, yielding a 2-way range ambiguity interval 
of approximately 150,000 km. 

The majority -vote combining logic in the old-r system was 
chosen [3] , [4] because it evenly (and optimally, for that strat- 
egy) distributed the power among all components for sequen- 
tial detection. All component delay measurements were thus 
made with times approximately proportional to the compo- 
nent periods. The clock-code (period 2) phase was acquired 
first, and then the other 75 code phase correlations were made 
sequentially. Two correlator channels were time-shared during 
range acquisition in a way that balanced the channel gains and 
removed any residual dc bias voltage in the baseband detection 
process. 

The reason for utilizing only a few correlators in ranging 
receivers until now has been that each correlator channel con- 
sisted of relatively expensive analog and unit-logic digital 
hardware. However, with the advent of reasonably high speed 
digital analog-to-digital devices and very large scale integrated 
(VLSI) digital devices, correlators may now be made at much 
more modest cost [7] . 


tern is now configured, about half of the range-measurement 
time is spent in correlating with the highest-frequency compo- 
nent (the clock), the other half being spent determining the 
range-cells of the lower-frequency components. The in system 
thus outperforms it by some 2.5 dB in signal-to-noise ratio. 

A companion paper by the author and J. R. Smith [7] dis- 
cusses the requirements and conceptual design of the correla- 
tor channels and VLSI devices required. 


II. The Transmitted Code 

The code component periods remain the same as in the 
earlier r system: 2, 7, 1 1, 15, 19, and 23. The combining logic 
for the transmitted code, however, is now taken to be 

x(c) = c l XOR AND (c) (1) 

in which c = (c x , . . . , c 6 ), c { is the clock (period 2) sequence, 
XOR is the exclusive-OR function, and AND ( ) is the logi- 
cal-AND of all 6 component sequences. The logical-AND, or 
“unanimous-vote” logic, is the limiting case of the majority- 
vote logic used previously. 

Since AND (c) contains only one “1” in its truth table of 
64 entries, the in-phase cross-correlation of x(c) with c x will 
be [3] 


_ (64-l)-l 
**1 64 UV/ 


( 2 ) 


The in-phase cross-correlation of x(c) with a code of the form 
(c l XOR c f -), for i = 2, . . . , 6, is only 

R x, ml ir =0031 (3) 


It is now economically feasible under VLSI technology to 
put the needed number of correlators into a receiver to build a 
detector for each of the components at each symbol-delay of 
the composite code. (It still may be impractical to build a full 
matched filter for the overall transmitted code, however.) It 
is therefore appropriate to reevaluate the relative merits of the 
two ranging methods under the removed constraint that pre- 
viously fixed the number of correlator channels. 


(These figures are only approximate, being influenced slightly 
by the sense of imbalances between 0’s and l’s in each of the 
pseudonoise sequences. These imbalances can be chosen to 
further optimize the reception, but this is left as an exercise 
for the implementer.) 

The cross-correlations as functions of ranging delay are 
shown in Fig. 1 . 


This article shows that a new-r composite-component up- 
link code (or “rr,” for brevity) which utilizes a new combining 
logic for the transmitter code and a 77-correlator receiver is 
again favorable in performance. In fact, it is shown that the 
vr method is only about 0.25 dB below the performance of a 
matched filter for the optimal transmitter code. As the ju sys- 


III. Performance 

The vr receiver achieves range measurement precision by 
clock-component correlation, just as did both predecessors. 
But since the clock component of the transmitted code con- 
tains 94 percent (i.e., 0.97 2 ) of the total ranging power, there 
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is only a 0.27 dB degradation in acquisition time from trans- 
mitting the clock component alone (such a code would not, 
however, remove the ambiguity of the range). 

The requirements for range accuracy demand that the stan- 
dard deviation of the range measurement due to noise be 
about 1/1000 of the “chip” time, t 0 , or symbol-rate of the 
clock. The relative clock variance o\ of A = r/t 0 is thus about 
10" 6 . In order to achieve this accuracy, the received signal- 
energy /noise-density ratio must accordingly be high [5] . 


The derivation of the maximum-likelihood detector is 
straightforward: We presume that we receive the transmitted 
binary ranging signal x(f) = x(c[f]), normalized here to unit 
power, immersed in wideband Gaussian noise n(t), as 

y(t) = otx(t - r) + n(t) (6) 

where ot = S 1 / 2 . The time delay t is to be estimated as that 
value ? maximizing the conditional probability (density) 
function 


ST 1 

67X10 (4) 


where S is the total signal power, T is the correlator integra- 
tion time, and N 0 is the received noise (single-sided) spectral 
density. 

Detection of the pseudonoise component phase is accom- 
plished by correlating the received signal with each of the 
separate phases of {c x XORc,). The power in each component 
is only about 0.001 = 0.31 2 of the total, so the component 
detection-energy/noise-density ratio is about 


S.T R 2 .ST 

~~j ~. — = x ‘ r = 0.001 X 6.7 X 10 4 = 67 (5) 


p{T\y(t),0<t<T} (7) 

Under the usual assumption that r is uniformly distributed 
over the unambiguous-range interval, the likelihood ratio, by 
Bayes’ rule, is 


, _ P{r I.K0} _ P{y(t) It} ^ 


( 8 ) 


where the interval (0, T) dependency has been suppressed for 
notational convenience. 

The probability of receiving y(t), given r, is the probability 
that the noise in n(t) = j'(r) - ooc(t - r). Because of the wide- 
band Gaussian character of the noise, the likelihood ratio 
becomes [6] 


The maximum required pseudonoise component-energy/ 
noise-density ratio for an error probability of 0.01 was about 
10 for the old-r system. Since detection of the pseudonoise 
range cell (see below) involves summing channel values for 
pairs of correlators, the additional noise may degrade the 
required S f T/N 0 to about 20 (a full analysis has not yet been 
made). The better-than-a-factor-of-three margin ensures, how- 
ever, that unerring full-range acquisition is almost certain. 




> 1 ( 9 ) 


where exp(x) is the exponential function e x 


IV. Maximum-Likelihood Receiver 

The maximum-likelihood estimator of the clock compo- 
nent phase is well known and will not be repeated here. The 
remaining pseudonoise code delays, however, are to be esti- 
mated from measurements made in parallel with the clock 
delay determination. This approach differs significantly from 
the old-r method, where the receiver pseudonoise sequences 
were acquired after the clock so that the receiver codes could 
be adjusted to then be in step with the received signal. Thus, 
whereas the old-r system enjoyed full component correlation 
in only one integration bin per component, the vt scheme 
must make do with partial component correlation in two adja- 
cent correlation channels for each component. 


By noting x 2 {t) - 1, canceling like terms in the numerator 
and denominator, and taking logarithms, we find that the con- 
dition on r is that 

t ~r 

y(t)x(t-'r)dt> I y(t)x(t-r)dt (10) 
\ Jo 

That is, r will be the maximum-likelihood estimator of r 
provided that it maximizes the correlation between the ob- 
served y(t) and the delayed transmitted code. However, a 
continuum of correlators is infeasible, so we must infer? from 
the finite number of measurements we do make. It has been 
shown [3] that the maximum-likelihood value can be inferred 
using correlations of the incoming signal and various delays of 
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the transmitted codes, in the form x y = c { and x t = c x XOR 
c t for i = 2, , 6. 



?.)* 


on 


out-of-phase correlation levels. Let z ( y denote the normalized 
sum of adjacent correlators at the /th position of the ith code 
component: 


, JVW 

Z H a. 


To decrease ranging inaccuracy caused by waveform distor- 
tion within the communication system, estimation of the 
clock component phase in the current (i system is performed 
by maximum-likelihood methods applied only to the funda- 
mental harmonic of the received clock component. This 
results in a modest, justifiable increase in required integration 
time. The vr method would presumably have the same require- 
ment for this clock estimation scheme. 


= R X i( T -f t 0 ) +R xi ( T -0'^)t o ) + n ii 


z.. + n.. 
i] ij 


(13) 


for an appropriately defined noise term The signal portion 
of z>j, denoted z { . in the equation above (see Fig. 2), 


Since almost all of the transmitted power is in the clock 
component of the code, the contribution of the 75 other cor- 
relators in improving the accuracy of the clock phase estimate 
will be negligible. Hence, the value of (r mod 2 t Q ) may be 
determined from the clock-channel correlators alone. (Com- 
bined, weighted estimation of the clock phase from all chan- 
nels can be done, however, with only slightly more complexity 
in the estimation program, if desired.) Since 94 percent of the 
transmitted power is in the clock component, and since 
maximum-likelihood estimation is performed on this compo- 
nent, the clock phase estimate is very nearly the same as the 
maximum-likelihood estimate of a pure clock signal. 

Thus, any method that with high likelihood selects the 
proper range-cell delays of the remaining components will 
measure the range within 0.27 dB of the performance of a 
maximum-likelihood device, insofar as ranging accuracy is 
concerned. 

The clock-channel measurement of (r mod 2t 0 ) is required 
to be very close to the actual value of (r mod 2t 0 ) for system 
accuracy. From this value, (f mod t 0 ) may be determined, as 
well as the ±1 sense of the in-step correlation (Fig. 1). There- 
fore, determination of (r mod Nt 0 ) additionally requires only 
the estimation of the integer values k { such that 

T-k.t 0 = (? mod t Q ) (12) 

for each of the remaining pseudonoise components (period 
N t t 0 ). We may estimate these values from the correlator out- 
puts of each pseudonoise component and combine them to 
form the overall range using the Chinese remainder theorem, 
just as did the previous r system. 

Only two of the correlator integration values/., /= 1, . . . , 
N t for the /th code component may derive from partial corre- 
lation with the true delay. The other I tj values correspond to 



If k is the correct range cell, i.e., kt Q < r < (k + l)r 0 , then 
the geometry of the correlation function (Fig. 1) leads to the 
estimator 


i 1 ik 


V*.-{ srfiwj^ (,5 > 


The right-hand side of this equation should be (tj mod f 0 ), 
as already estimated accurately by clock -channel computa- 
tions, within the expected noise deviation. 

If there were no noise, the value of j yielding the maximum 
z t j would be the correct range cell. But with noise, the cells on 
either side must also be scrutinized, as the spillover correlation 
and noise may cause us otherwise to choose the index of the 
wrong range cell. 

Hence, let k i be thay‘-index for which z (j is maximum, and 
let k ( be 'k i ~ 1 , k h or k f + 1 , whichever minimizes the differ- 
ence between and ? l . This k ( is a high-likelihood range- 
cell estimate for the zth component because the maximum- 
likelihood range cell is certainly one of the three candidates, 
and the comparison of the symbol-fraction range offset, calcu- 
lated as above, against the accurate determination of the clock 
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symbol fraction, rules out the other two candidates with high 
probability. 

Analysis indicates that the probability of making the wrong 
choice would be very remote. As may be seen in the geometry 
of the cross-correlation function, shown in Fig. 1, the wrong 
pair of correlator values inserted into the symbol fraction 
formula above produces a significantly different estimate for 
r f . than for Tj. 

V. Range Calculation Algorithm 

The procedure used by the i n receiver to reconstruct the 
range is therefore basically the same as in the old-r system, 


except for the way the component range cell determinations 
are made. Having read the 77 correlation values all at once 
after an integration interval T, the normal calculations on the 
two clock values determine the clock component range delay 
(mod 2t 0 ) with high accuracy. The remaining 75 correlation 
values are grouped by component, and then a value k t is 
chosen for each component. This k i is the range-cell index / 
that maximizes the sum of adjacent correlation values I t j 
+ Each of the three range-cell indices k i - 1, k h and 

k f + 1 in turn is hypothesized to be the correct range-cell in- 
dex. The appropriate correlator values for each of these indices 
are then used to compute a symbol-fraction offset using the 
estimator formula above. The candidate that comes closest to 
the clock-channel measurement wins. 
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Fig. 1. Composite-code cross-correlation function R x j (x) 



RANGE DELAY 

Fig. 2. Sum of adjacent correlator outputs versus range delay 
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